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INTRODUCTION 

The  ei-a  subsequent  to  the  last  World  War  has  seen  tremendous 
technological  strides  in  the  fields  of  nianned  and  unmanned  flights. 
Manned  flights  of  aerodynamic  and  hybrid  vehicles  have  extended  to 
speeds  of  several  times  that  of  sound  and  to  altitudes  of  tens  of  miles. 
Control  of  these  vehicles  has  long  since  passed  the  stage  of  direct  con- 
trol and  fly  by  feel.     Very  sopisticated  control  systems  are  used  to  aid 
the  pilot  in  manned  flight  and  to  completely  control  the  vehicle  in  the 
case  of  unmanned  flights. 

One  of  the  important  subsystems  in  these  sophisticated  flight 
control  systems  is  the  automatic  pilot  or  autopilot.     It  is  the  function 
of  the  autopilot  to  carry  out  precisely  the  maneuvers  called  for  by  the 
guidance  system.      Coupled  together,    the  autopilot  and  guidance  system 
direct  the  vehicle  along  the  trajectory  required  to  intercept  the  target. 
Accuracy,    speed,    and  stability  are  three  prime  requirements  for  a 
good  autopilot. 

In  the  missile  field,    vehicles  of  extremely  high  speed  and  altitude 
capabilities  have  been  developed.      These  vehicles  require  control  sys- 
tems which  are  extremely  fast  and  which  have  very  broad  dynamic 
ranges.      To  meet  this  severe  demand  on  the  control  system  and  in 
particular  on  the  autopilot  the  use  of  intentionally  introduced  nonlinear 
elements  is  often  considered. 

A  typical  nonlinear  element  is  one  which  has  only  two  possible  states, 
fully  on  in  either  of  two  directions.     Such  devices  usually  simplify  the 


system  hardware  but  almost  always  complicate  the  analysis.     Super- 
position is  no  longer  permissible  and  analysis  for  one  type  of  signal 
does  not  always  yield  useful  inforination  about  other  types.      The 
merits  of  such  systems,    however,    more  than  offset  the  difficulty  of 
analysis.  -^  ' 

The  analysis  of  a  simple  nonlinear  control  system  for  an  aero- 
dynamic vehicle  is  undertaken  in  this  paper.     The  prinnary  objectives 
of  the  work  to  follow  are  to  show  a  means  of  analyzing  a  system  of  this 
type  and  to  indicate  a  design  procedure  for  a  fixed  parameter  system. 
For  a  simple  system  amenable  to  a  direct  literal  analysis,    a  complete 
design  is  undertaken.     It  is  shown  that  air  data  scheduling  is  not 
required  to  achieve  a  stable  system,  with  adequate  response  throughout 
the  flight  regime  if  a  two  state  nonlinearity  is  used. 

For  any  aerodynamic  configuration  the  effectiveness  of  the  con- 
trol surfaces  in  changing  the  flight  path  is  a  function  of  the  Mach  number 
and  the  local  static  pressure.     Hence,    the  control  system  must  cope  with 
a  vehicle  whose  equations  of  motion  are  strong  functions  of  the  local 
dynamic  pressure.     Typical  flights  may  encompass  a  range  of  dynamic 
pressure  of  ten  to  a  hundred  depending  on  the  vehicle  and  the  trajectory. 
This  fact  alone  precludes  the  use  of  a  fixed  parameter  linear  system. 
As  a  consequence,    early  linear  control  systems  compensated  for  dynamic 
pressure  variations  by  elaborate  use  of  air  data  measurements.     Mach 
n\imber  and  local  static  pressure  were  continuously  measured  (usually 
by  measuring  ram  pressure)  and  these  data  used  to  vary  the  control 
system  parameters. 


To  design  a  system  of  the  above  type,    one  must  have  accurate 
knowledge  of  the  airframe  transfer  function  for  all  dynamic  pressures 
and  all  trim  conditions.     Furthermore,    the  hardware  necessary  to 
accomplish  the  gain  scheduling  is  quite  connplicated  and  often  un- 
reliable.     Use  of  nonlinear  control  as  suggested  in  this  paper,    obviates 
the  need  for  air  data  scheduling  equipmSht  and  at  the  same  timie  yields 
an  autopilot  which  has  a  response  superior  to  a  linear  system. 

SYSTEM   DESCRIPTION 

The  functional  block  diagram  of  a  typical  guided  missile  control 
system  is  shown  in  Fig.    1.     In  any  such  system  the  target  and  missile 
positions  in  space  are  compared  and  the  resultant  guidance  error  signal 
is  fed  to  the  missile  control  system  which  in  turn  acts  to  null  the  error 
signal.      Clearly,    the  particular  method  of  guidance  error  detection  is 
determined  by  the  type  of  guidance  system  used.     For  example,    a  beam- 
riding  systenn  eniploys  a  separate  guidance  and  tracking  radar  whereas 
a  homing  system  employs  a  missile  bound  seeker  unit.      The  type  of 
guidance  system  used  will  also  determine  the  basic  nature  of  the  guidance 
error  signal. 

The  guidance  computer  is  used  to  provide  the  appropriate  error 
signal  gain  and  shaping  to  convert  the  error  signal  to  a  suitable  missile 
lateral  acceleration  command  signal;     thus  its  form  is  determined  by 
the  type  of  guidance.      The  missile  autopilot  in  turn  accepts  the  acceler- 
ation command  signal  and  activates  the  control  surfaces  as  required  to 
provide  this  desired  lateral  acceleration.     The  component  of  this  lateral 


acceleration,   which  is  normal  to  the  flight  path  acting  in  conjunction  with 
the  missile  velocity,    produces  a  change  in  flight  path  angle  so  as  to  null 
the  error  between  the  desired  missile  position  in  space  (  the  target 
position)  and  the  actual  position. 

That  portion  of  the  missile  system  which  constitutes  the  autopilot 
is  enclosed  in  dashed  lines  in  Fig.    1.     It  is  this  portion  of  the  system 
which  will  be  of  concern  for  the  remainder  of  this  paper.      The  airfranae 
considered  in  this  study  uses  small  control  surfaces  located  remotely 
from  the  missile  center  of  gravity  to  achieve  large  turning  naoments. 
Lift  or  lateral  acceleration  forces  are  derived  primarily  from  the  missile 
body  and  thus  are  due  primarily  to  angle  of  attack.     Normal  forces  result- 
ing from  control  surface  deflection  are  small  relative  to  the  body  forces 
and  are  thus  neglected. 

The  airborne  portion  of  a  missile  control  system  must  maintain 
control  of  the  three  airframe  coordinates;    pitch,    yaw,    and  roll. 
Generally,    the  roll  system  functions  only  to  stabilize  the  roll  angle.     A 
roll  free  gyro  and  a  roll  rate  gyro  are  used  in  conjunction  with  appro- 
priate feedback  networks  to  stabilize  the  roll  angle  (i.e.  ,    the  roll  system 
maintains  a  prescribed  roll  orientation  throughout  the  flight).     For  pur- 
poses of  this  study,    it  is  assumed  that  separate  moment  or  force  devices 
are  provided  for  the  roll  system.     This  essentially  divorces  the  roll 
and  steering  systems  so  that  they  may  be  treated  as  distinct  systems. 

The  pitch  and  yaw  planes  are  essentially  identical  and  independent 
in  the  absence  of  roll-yaw-pitch  (RYP)   coupling  through  the  aerodynamics 
and  thus  treatment  of  one  plane  suffices  for  both.     Maneuvers  in  either 


plane  or  in  any  connbined  plane  are  permissible  with  the  latter  being 
a  linear  combination  of  the  single  plane  motions.     Separate  sensing 
instr undents  are  provided  for  each  plane. 

Standard  pitch  plane  nomenclature  is  shown  in  Fig.    2.     In  the 
discussion  to  follow,    the  ternns  pitch,    yaw,    and  steering  refer  to  lateral 
motions  or  planes  of  motion  and  are  used  synonomously.     In  Fig.    1, 
the  autopilot  is  assumed  to  be  an  "acceleration  feedback"  (AFB)    system 
wherein  the  command  acceleration  is  compared  with  the  actual  missile 
acceleration  as  measured  by  an  accelerometer .     The  acceleration  error 
thus  obtained  is  used  to  activate  the  control  surfaces  to  null  the  error. 
Missile  attitude  rate  feedback  is  required  to  provide  sufficient  autopilot 
loop  damping.     Acceleration  rate,    control  surface  position,    and  control 
surface  rate  are  also  often  used  as  feedback  quantities  in  system 
synthesis. 

The  above  description  is  of  a  typical  high  speed  aerodynamic 
vehicle  and  control  system.     The  unique  feature  of  the  system  con- 
sidered in  this  paper  is  the  control  surface  actuator.     Rather  than  con- 
trolling surface  position  or  surface  rate  in  direct  proportion  to  an  actu- 
ation signal,    this  system  is  controlled  by  applying  full  tail  (control  sur- 
face) rate  at  all  times.     This  fact  is  expressed  by  the  control  equations: 

6    =  +B       m  >  0  (la) 

6   =  -B       m  <  0  (lb) 

This  method  of  control  gives  rise  to  some   special  types  of  behavior 
which  are  discussed  in  the  next  section. 


STABILITY   ANALYSIS 
Qualitative  Discussion 

Figure  3  shows  a  simple  diagram  of  this  nonlinear  system  with  a 
symbolic  representation  of  the  actuator.     In  this  system,    G(s)  and  H(s) 
are  such  that  the  autopilot  loop  is  conditionally  stable  (i.e.  ,    a  linear 
actuator  coupled  with    G(s)  and  H(s)  would  yield  a  conditionally  stable 
system).     Closing  a  conditionally  stable  loop  with  a  bi- stable  element 
such  as  in  this  case  gives  rise  to  at  least  two  modes  of  oscillatory 
operation.      These  are  a  low  frequency  unstable  oscillation  and  a  higher 
frequency  stable  oscillation.     Under  most  situations,    the  system  will 
converge  to  the  stable  oscillating  point. 

For  manned  aircraft,    an  oscillating  autopilot  is  generally  regarded 
as  unacceptable  for  several  reasons.     Jitter  of  the  control  stick  result- 
ing from  the  oscillation  is  objectionable  to  the  pilot  and  wear  on  the  air- 
craft systems  due  to    continual  chatter  is  undesirable.     Missiles,    how- 
ever,   have  a  very  short  mission  time  and  are  nearly  always  one  shot 
applications  so  wear  and  fatigue  are  not  problems. 

The  higher  frequency  stable  limit  cycle  oscillation  is  the  mode  of 
operation  desired  in  this  control  system..     A  linait  cycling  state  results 
in  an  approximately  linear  characteristic  for  the  actuator  for  signals 
which  do  not  exceed  the  limit  cycle  amplitude.     For  snaall  signals  then 
the  loop  may  be  analyzed  using  linear  theory  and  a  pseudo-linear  gain 
term  for  the  bi- stable  element.      The  theory  of  such  analyses  is  well 
documented  by  Gorozdos  (z)  and  is  not  covered  in  any  detail  in  this  report. 


It  should  be  noted  that  for  very  small  amplitude  oscillations  a  linear 
analysis  is  of  little  use  since  the  oscillation  is  suppressed  throughout 
the  major  portion  of  most  transients. 

Large  inputs  or  disturbances  drive  the  system  operating  point 
toward  the  lower  frequency  oscillation.     If  this  point  is  ever  reached 
the  system  oscillates  in  a  diverging  manner  with  frequency  approaching 
zero   and  amplitude  approaching  infinity.     Hence,    this  lower  critical 
frequency  is  a  catastrophic  operating  point  which  the  system  must  never 
be  allowed  to  reach.      Unless  suitable  design  precautions  are  taken,    this 
factor  represents  a  limit  on  the  permissible  command  level. 

For  intermediate  level  inputs  which  are  of  most  interest  from  a 
response  point  of  view,    the  actuator  is  not  even  approxinriately  linear  and 
n'lust  be  treated  as  a  bang-bang  device.     When  nonlinear  control  of  this 
type  is  used,    there  exists  what  is  called  an  optimiim  bang-bang  solution 
or  response.     If  such  a  response  can  be  achieved,    the  error  in  response 
to  a  step  input  is  reduced  to  zero  in  the  minimum  possible  time. 
Methods  appropriate  to  this  type  of  solution  are  applied  in  a  later  section 
to  yield  the  optimum  form  of  the  response  and  the  appropriate  data  for 
synthesis  of  an  optimum  systenn.      '. 

The  small  signal  and  stability  problems  are  handled  by  developing 
a  set  of  stability  surfaces  from  describing  function  methods.     These  sur- 
faces delineate,    in  a  conveniently  usable  fashion,    all  the  possible  points 
of  instability  in  a  three  dimensional  parameter  space.     By  combining 
surfaces  for  various  operating  regimes,    a  set  of  fixed  parameters  most 
suitable  to  the  entire  flight  regime  may  be  selected.      These  surfaces  also 


indicate,    in  a  very  vivid  way,    the  influence  of  such  factors  as  airframe 
stability  and  actuator  paranneters  on  overall  system  stability. 

Stability  Criterion  and  Characteristic  Equation 

Figure  4  shows  a  block  diagram  of  the  specific  control  system  which 
is  analyzed.     In  the  absence  of  inputs  or  initial  conditions,    a  system  with 
a  nonlinearity  of  the  type  shown  in  Fig.    4  may  exhibit  any  one  of  several 
modes  of  behavior.     To  determine  the  nature  of  the  oscillating  modes  of 
operation,    it  is  necessary  to  apply  a  suitable  stability  criterion.     For 
stability  analyses  of  systems  with  simple  nonlinear  elements,    like  the 
switch  under  consideration,    a  describing  function  is  a  satisfactory  ap- 
proximation to  the  true  switch  characteristic.     A  switch  with  output  limit 
B  and  input  sine  wave  of  amplitude    a    has  the  describing  function: 

■■■-.    '    ■-■      ff     ■ 
N  =    4B/Tra  -      •"'  (Z) 

Notice  that  the  effective  gain  of  the  element  decreases  as  the 
input  level  increases.     This  is  characteristic  of  saturating  ele- 
ments,   and  is  of  monumental  consequence  with  regard  to  condi- 
tionally stable  systems  as  the  discussion  to  follow  will  show. 

With  the  nonlinear  element  represented  by  its  describing 
function  N,    the  simplest  form  of  the  autopilot  loop  is  sh6wn  in  the 
following  sketch.      The  approximate  characteristic  equation  for  this 
loop  is: 

1  +  NG(s)   H(s)  =  0  (3) 


n. 


This  loop  will  oscillate  at  the  angular  frequency    oj        provided: 


G(jco^)   H(jo3j   =   -1/N 


(4a) 


Equation  4a    may  be  equivalently  written  as: 


(4c) 


iG(joj^)  H(jcoJ  1   =  1/N  (4b} 

/G(jco^)   H(jco^)     =  ±180° 

If  G(s)  and  H(s)  are  specified,    then  any  of  the  graphical  methods  of 
analysis  (Bode  plots,    Nyquist  plots,    etc.)  will  suffice  to  determine 
the  possible  frequencies  and  amplitudes  of  oscillation.      With  these 
methods,    however,    much  effort  is  required  to  find  a  single  solution 
of  Eq.   (3).       One  has  only  experience  and  past  computations  as  a  guide 
to  selecting    G{s)  and  H(s).      Considerable  trial  and  error  effort  is  in- 
volved in  arriving  at  a  design  which  even  then  may  not  be  optimum. 
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An  alternate  procedure  is  to  solve  Eq.    (3)   explicitly  for  the  points 
of  intersection  of  l/N  and  G(s)  H(  s)   rather  than  plotting  the  two  and 
finding  the  points  graphically.      This  is  accomplished  by  assuming  an 
oscillation  at  the  frequency    w        and  then  solving  for  the  free  para- 
meters in    G{s)  and  H(s}    necessary  to  insure  this  oscillation.     These 
computations  are  carried  out  very  quickly  by  machine  methods,    although 
hand  methods  would  yield  useful  data  in  a  reasonable  time. 

For  a  steady  oscillation  to  exist  at  least  one  pair  of  roots  (a  con- 
jugate pair)  of  the  characteristic  equation  must  lie  on  the  imaginary 
axis.     Hence,    the  necessary  condition  for  an  oscillation  at  the  angular 
frequency    oj      is:  .  - 


o 


1+  NG(jco^)   H(jco^)  =  0       .  (5) 

Separating  this  equation  into  real  and  imaginary  parts  one  gets: 

U(N,w^)  +  jV(N,  oj^)  =  0  (6) 

The  real  and  imaginary  parts  must  be  zero  simultaneously  and 
hence,    Eq.   (6)  provides  the  two  equations: 

U(N,   00^)  =  0  (7a) 

V(N,   w^)  =  0  (7b) 

These  equations  could  be  solved  for  the  frequency  of  oscillation 

w^    and  the  corresponding    N    for  a  particular  set  of  loop  parameters. 

In  actual  practice,    however,    the  equations  are  transcendental,    and  it 
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is  niore  convenient  to  write  them  in  the  form: 


U'(l/NK^,    ^x^^n^  "  °  ■  ^^^^ 

V'(1/NK       K    /K    )  =  0       ■  (8b) 


n        X       n 


By  specifying    w     ,     these  equations  may  be  solved  for  the  pair 

(l/NK    ,    K    /K     ).       In  this  form    K       denotes  any  one  of  the  feed- 
n        X       n  X  ' 

back  gains,    and  all  parameters  are  normalized  with  respect  to 
the  scale  factor    K    .       In  solving  for  a  specific  pair  of  parameters, 
all  others  are  held  constant.       A  single  solution  of  Eqs.   (8)  yields 
a  point  (l/NK^,    K^  /K^,  K^   /K^,    K^.  /K^  )  in  four  space  with  run- 
ning parameter    oj     .       Later  argmnents  will  show  that  relative 
values  of    K.     and    K-       are  constrained  by  the  flight  regime.      Con- 
sequently,   it  is  simpler  and  more  convenient  to  define  the  ratio 

^n  /^e     =     ^r  ^  -  -  (9) 

and  then  to  use  Eqs.   (8)   to  develop  a  surface  in  the  three  space 

with  coordinates    (l/NK       K^   /K    ,    Kv  /K     ).       This  surface  contains 

n         o         n        0         n 

w       as  a  parameter  and  is  developed  for  a  fixed  value  of    K   . 

Derivation  of  the  Stability  Equations 

To  find  explicit  stability  surfaces  for  the  autopilot,    it  is  necessary 
to  derive  Eqs.  (8).       Transfer  functions  for  the  airframe  under  various 
flight  conditions  are  presented  in  Appendix  A  of  this  paper.     The  accel- 
erometer  and  gyro  are  assumed  to  be  characterized  by  half  critically 
damped  quadratic  terms.      The  control  surface  actuator  is  an  open  loop 
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device  with  a  pure  transportation  delay  of    t    seconds.       The  block 
diagram  of  the  composite  system  showing  all  of  these  features  is 
shown  in  Fig.    5. 

Figure  5  is  easily  reduced  to  the  single  loop  system  of  Fig.    6. 
Distinct  parameters  for  the  accelerometer  and  gyro  are  used  in  the 
derivation  of  the  stability  equations  to  permit  study  of  the  influence 
of  various  instruments  on  stability.     In  general,    these  instruments 
have  slightly  different  parameters,    although  they  have  been  assumed 
to  be  alike  for  the  numerical  results  presented  in  this  paper. 

Substituting  directly  from  Fig.    6  into  Eq.    (3)  and  simplifying 
gives: 

as''    +  {^  +  a  AK)s'^+  {7+  ^  AK  -  Ca  )  s^  +  ^  +  AK7  -  C/3)  s"^ 

+  (l  +  AKs^   =  CI)  sH  (AK  -   C^)s2  -   Cs  +  Ne'^^  JK-    DT|    s^ 

+(Kg  ETi  -f  K^Dt|+  2K^DCgTg)s^+  (K^D+  2C  ^T^K^' E -i-  2C  ^T^K^D 

+(K^Afg+  Tl  {K^E    K-^  KA}  )  s3+  (Kg  E  +  K^D  +  K^AT|    K'  KAT^^ 

+  '^a^a^^e^+  ^e'^^>+  'SVn^   s^+(  {K- E+ Kg  KA+K^A 

The    D    term  defining  the  normal  acceleration  per  unit  of  control 
surface  deflection  is  carried  through  the  literal  analysis  for  gener- 
ality but  is  dropped  in  the  numerical  analysis  because  of  its  small 

magnitude.     Equation  (lO)   is  to  be  solved  for     s  =     iw        thus  allowing 

o  ° 

the  exponential  to  be  expanded  according  to  the  relation: 

-JGO       T 

e        o      =cosw     T-     j  sin  w    t  -  / 1  ij 
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Substituting     s  =  joj        into  Eq.   (lO),    using  Eq.   ( 1 1)  ,    and  separating 

into  i-eal  and  innaginary  parts  gives  the  following  pair  of  stability 

equations:  .         •  .. 

-(jw^  +  (7  +  /3  AK-   Ca)oj5  -  (l  +  AK^    _   C7)w^  -   Ceo 
o       ^  ^ '    o       ^ o o 

K    N 
n 


r    0 


K^      [^  ^'g    o 


|jl(T^oo^-1)u>       cosco      t  -   2i     Too^sinw      T 


g    g    o 


K-. 


+  -^ —       (KA-7Toj^)a)      cosw     T-((i)co^-9)w^sinw     t 
K'  o'o  o^^o'o  o 

n      U-  — 


K- 


K 


n 


1 


-(cbco^    -9)w      cosoj     T+(-iTw^    -KA)    sin  CO     T 
^^     o        '      o  o       ^        o  o 


r 

-    (Zui^     T     CO       cosco     T     fx{T^co^    -    l)    sin  co     T 
[ggo  o'^'go  o_ 


=  0 


(12a) 
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For  any  parameter  set  and  frequency    co       which  satisfies  Eqs.   (iZ), 
the  system  will  oscillate  at  the  angular  frequency    co    .     Further 


(12b) 
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tests  are  necessary  to  determine  if  the  oscillation  is  stable  or  unstable. 
This  topic  is  reserved  for  discussion  at  a  later  point  in  this  paper. 

Solution  of  Stability  Eqiiations 

Examination  of  Eqs.   ( I  2)   shows  that  upon  substitution  of  a  numer- 
ical value  for    w     ,     the  equations  become  algebraic.     Note,    however, 
that  there  are  only  two  equations  in  three  unknowns  even  if    K       is  held 
constant.     A  solution  is  obtained,    however,    by  selecting  one  parameter 
(usually    K-    or  Kv     )  and  solving  for  values  of  the  remaining  pair  for 
all  w       which  have  naeaning.     Obviously,    only  real  positive  solutions 
are  of  interest.     Furthermore,    one  must  exercise  care  in  selecting 
values  of  the  parameter  to  be  held  constant  to  insure  that  a  solution 
for  the  other  pair  of  parameters  exists.     For  machine  computations, 
Eqs.   (12)  are  written  in  the  general  forms 

(I3a) 
A2(<.^)  /K^+  K^K-  B2(c^J  /K^  +  K-    C,  (co^)  /K^  +  K^-  Dz{u^^)  /K^+  £2(0;^)  =  0 

:      .       -        ,  .      ■  (13b) 

where  the    Ai(a)^) ,    Bi(co^)  ,    A2(co  ) ,    etc.,    are  algebraic  or  trans- 
cendental functions  of    co  .     For  specific  values  of    co  ,    these  are 

o  ^  o 

ordinary  algebraic  equations  which  may  be  solved  using  Krameres  rule 

or  any  other  nnethod  for  simultaneous  equations. 

Every  solution  of  Eqs.   (l3)  yields  a  point  in  three  space  with 

coordinates  (l/KN,    K-     /K    ,     K'-    /K    )   and  with  characteristic 
n  o  n  u         n 
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frequency    to  .       Repeated  solution  of  Eqs.   ( 1  3)  will  thus  yield  a  sur- 
face in  three  space  with  frequency  as  a  parameter.     This  surface 
defines  all  possible  states  of  oscillation  for  the  given  system  and  thus 
is  a  very  useful  tool  in  selecting  an  operating  point  and  in  analyzing 
the  system  for  small  inputs.      It  is  now  necessary  to  determine  the 
exact  meaning  of  these  surfaces  in  terms  of  system  operation. 

The  basic  plant  and  feedback  functions  of  the  system  under  con- 
sideration would  yield  a  conditionally  stable  system  if  the  nonlinear 
element  were  replaced  with  a  linear  one.     Hence,    a  stable  high  fre- 
quency oscillation  and  an  unstable  low  frequency  oscillation  are  pos- 
sible.    To  investigate  the  stability  of  these  oscillations,    one  may  con- 
sider the  nonlinear  element  as  a  linear  device  whose  gain  decreases 
with  increasing  input  level  (see  Eq.   (Z)).     By  observing  the  motion  of 
the  roots  of  the  characteristic  equation  on  the  s-plane  as  the  oscilla- 
tion annplitude  changes,    one  may  determine  the  stability  of  a  particular 
oscillating  point.     (For  purposes  of  this  analysis  a  stable  oscillation 
is  one  which  has  a  constant  amplitude  and  frequency  and  which,    when 
perturbed,    returns  to  this  point.     An  unstable  oscillation  is  one  which 
violates  any  or  all  of  these  conditions.)     Consider  the  root  locus  sketch 


of  a  conditionally  stable  shown  below 
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At  the  higher  oscillation  frequency    u:^,     a  positive  perturbation  in 
amplitude  (which  decreases  the  switch  gain),    causes  the  operating 
point  to  shift  into  the  left  half  plane.     Here  the  system  is  stable  and 
the  oscillation  tends  to  die  out.      This  decay  in  amplitude  increases 
the  gain  of  the  nonlinear  element  and  the  operating  point  converges 
on    0^    as  a  stable  operating  point  (an  analogous  argument  to  that  given 
above  will  show  that  a  perturbation  into  the  right  half  plane  at    cj^    will 
also  converge  to    ci^  . 

Very  large  signals  appearing  at  the  switch  input  cause  the  oper- 
ating point  to  shift  from    c*^    toward  the  origin  along  the  locus.     Assume 
that  the  switch  gain  is  driven  sufficiently  low  to  reach  the  point    cjj. 
At  this  point  a  slight  increase  in  the  switch  input  (causing  a  decrease 
in  gain)   shifts  the  operating  point  into  the  right  half  plane.      With  roots 
thus  in  the  unstable  part  of  the  s-plane,    the  amplitude  of  oscillation 
tends  to  grow.     This  further  decreases  the  gain  causing  the  systemi 
to  diverge  with    w      approaching  zero,    and  the  amplitude  of  oscilla- 
tion approaching  infinity. 

A  second  point  should  be  noted  here.     As  the  gain  changes  with 
varying  signal  level,    the  roots  of  the  characteristic  equation  move 
along  the  root  locus.       One  would  thus  expect  a  response  which  is   sensi- 
tive to  input  amplitude.      This  is  indeed  the  case  and  this  point  will  be 
explored  later  in  this  paper.     Furthermore,    from  the  above  discussion, 
one  must  beware  of  inputs  which  exceed  the  unsafe  boundary  at    coj.       The 
stability  surfaces  which  can  be  developed  from  Eqs.   ( 1  3) ,    delineate  the 
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possible  oscillation  frequencies  (both  the  stable  and  unstable  ones) 
and  also  show  in  a  qualitative  sense  the  permissible  disturbance  level 
for  stability. 

Interpretation  of  Stability  Surfaces 

Equations  ( 1  3)  are  solved  for  all  values  of    w      which  yield  real, 
positive  sets  of  parameters.       These  data  could  be  presented  as  a  three 
dimensional  plot  or  surface  with  coordinates  (l/K   N,    Ki  /K^,  K-^  /K^). 
■With  the  range  of  parameters  involved,    however,    it  is  mxuch  more  con- 
venient to  consider  planes  taken  from  the  surface  at  constant  values 
of    Kj.      or    K-A     as  is  illustrated  in  the  sketch  below. 


K-^- 


K 


n. 


K*    =  constant 


Kx    =  constant 


A  typical  intersection  of  a  plane  of  constant    K-      and  the  stability- 
surface  is  sketched  below.      The  direction  of  increasing  oscillation  fre- 
quency is  indicated  by  the  arrow  on  the  sketch.      This  result  may  be 
interpreted  by  referring  to  the  root  locus  sketch  for  a  conditionally 
stable  system  also  shown  on  this  page.     It  is  true  thai  the  root  locus 
method  m.ay  not  be  applied  directly  to  a  system  with  a  transportation 
delay.       ,,     ■  .  .;     '>  ' 
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However,    the  basic  form  of  the  stability  surfaces  is  not  altered,    since 
the  principle  influence  of  the  delay  is  to  change  slightly  the  frequency 
associated  with  a  particular  parameter  point.     Interpretation  of  the 
stability  boundaries  (  the  two  dimensional  planes  taken  from  the  surface) 
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in  terms  of  the  root  locus  is  thus  valid  in  a  qualitative  sense  although 
it  cannot  be  used  to  obtain  quantitative  data  unless  the  delay  is  zero. 

Returning  to  the  sketches  above,    operation  of  the  system  is 
restricted  to  the  path    A-B    on  the  stability  plane  and  also  on  the  root 
locus  plane.     The  nonlinear  element  may  be  treated  as  a  device  whose 
gain  varies  inversely  v/ith  the  input  amplitude  (see  Eq.   [Z]  ).      It  has 
previously   been    shown  that  a  temporary  operating  point  which  lies 
anywhere  between    A    and    A'    will  eventually  become  a  stable  oper- 
ating point  at    0    with  frequency    cj2.       However,    if  the  systenn  is  ever 
driven  to  points  in  the  region    A'  -  B,      the  system  will  diverge. 

Referring  to  the  sketches,    one  sees  that  decreasing    K?.      causes  more 

o 

and  more  of  the  root  locus  to  lie  in  the  right  half  plane.       When    K- 

8 

corresponds  to  point    C,      the  root  locus  is  tangent  to  the  imaginary 
axis  and  hence,    lies  entirely  in  the  right  half  plane.       With  these 
facts  in  mind,    the  following  interpretations  of  the  areas  shown  on  the 
stability  boundary  sketch  are  obvious: 

1.       Area    D    is  an  unstable  region  and  miust  be  ex- 
cluded as  an  operating  region.     This  area  may 
be  entered  by  setting    Ky      too  small  (Kv  <    C) 
or  by  subjecting  the  system  to  large  inputs  or 
disturbances  thus  traversing  the  boundary  along 
A-B. 
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2.  Area    E    is  a  safe  area  with  some  of  the  instant- 
aneous roots  being  in  the  left  half  plane.     Tempor- 
ary operating  points  in  this  region  will  converge 
eventually  to    G02    along    O-A'. 

3.  Area    F    is  also  a  safe  area,    but  the  instantan- 
eous roots  lie  in  the  right  half  plane.     Temporary 
points  in  this  area  also  converge  to    GJ2    but  here 
along  A-0 . 

4.  For    K"'      corresponding  to  point  A,    the  steady 
state  operating  point  is     0    with  oscillation  fre- 
quency   0J2* 

To  generate  a  boundary  of  the  type  discussed  above,    one  must,    of  course, 

select  a  value  of    Kt,     which  yields  real,   positive    K--  .       A  similar  state- 

ment  applies  to  boundaries  of    K-  . 

Planes  of  constant    K"     are  of  the  general  nature  sketched  below. 
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Using  arguments  similar  to  those  above,    one  finds  that  the  following 
interpretations  apply  to  the    K^     boundaries. 

1.  Values  of    K-      larger  than    D    result  in  the  entire 
root  locus  being  in  the  right  half  plane.     Hence, 
this  region  must  be  excluded  as  an  operating  region. 

2.  Any  value  of    K-      between    A    and    D    is  acceptable 
■with  the  lower  boundary  being  stable  and  the  upper 
boundary  unstable. 

3.  For  values  of    K-      less  than    A,     the  systenn,    in  an 
undisturbed  state,    will  exhibit  a  stable  high  fre- 
quency oscillation  at    coj.      This  linearizes  the  loop 
for  small  signals  and  hence  yields  a  conditionally 
stable  linear  loop.     The  low  frequency  gain  margin 
for  operation  in  this  region  is  negative,    and  the 
system  diverges.     This  divergence  is  characterized 
by  a  constant  amplitude  oscillation  at    co^    with  a 
superimposed  low  frequency  oscillation  with  ampli- 
tude approaching  infinity  and  frequency  approaching 
zero. 

Planes  of  constant     l/N    do  not  contribute  any  additional  information. 
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Niimerical  Resvilts 

Two  flight  regimes,    representing  a  typical  range  on   q*  ,      have 
been  selected  for  study.     In  addition,    an  extreme  range  of  airframe 
stability,    associated  with  various  center  of  gravity  locations,    is  con- 
sidered for  the  high    q    flight  condition.       These  combinations  are 
sufficient  to  illustrate  the  use  of  the  method  of  analysis  and  the  range 
of  parameters  to  be  expected.       Table  1-A  lists  aerodynamic  para- 
meters for  the  conditions  considered. 


Table  1-A.     Aerodynamic  parameters  used  in 
calculating  stability  data. 
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1.00 

0 

0 

229.0 

1.  1 

1 

1.00 

0 

-100  • 

229.0 

1.  1 

2 

1.00 

0 

100 

229.0 

1.  1 

3 

1.00 

0 

0 

229.0 

0 

4 

.  116 

0 

0 

26.0 

1.7 

5 

.  116 

0 

0 

26.0 

0 

6 

Table  1-B  lists,    for  the  aerodynamic  parameters  listed  in  Table 
1-A,    the  feedback  parameters  for  which  stability  data  have  been  ob- 
tained.    The  resulting  data  are  presented  in  the  figures  listed  in    the 
table  and  are  discussed  below. 


*     See  page  43  for  nomenclature. 
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Table  1-B.     Feedback  and  actuator  parameters 
used  in  calculating  stability  data. 


Flight  K-  K- 

condition 


e  e 


T 

Fig 

No. 

0.01 

7 

0.01 

8 

0.01 

9 

0.01 

10 

0.01 

11 

0.01 

12 

0.01 

13 

0.01 

14 

1  0.1      0.03      0.01 

5  0.1      0.03      0.01       ■" 

1  .  0.03      0.01         0.00  3 

5  0.01      0.003      0.001 

1,2,3  0.03 

1,  2,  3  0.003 

1  0.03 

1  0.003 


Figures  7  and  8  show    Kv     boundaries  for  high    q    and  low    q 

conditions,    respectively.       Similarly,      Figs.   9  and  10  are    K;      bound- 

0 

aries  for  high  and  low    q    conditions.       The  boundaries  associated  with 
the  different    q    levels  are  not  radically  different,    however,    it  should 
be  noted  that  the  high    q    condition  generally  dictates  the  values  of    K- 
and    Kg  .     The  actual  values  of    K-     and  Kr-      are  selected  based  on 
relative  stability,    small  signal  response,    and  to  some  extent,    the  time 
domain  analysis  presented  later  in  this  paper. 

Data  for  flight  conditions  1,    2,    and  3  are  presented  in  Figs.    H 
and  12.       These  data  show  that  airframe  stability  influences  the  upper 
portion  of  the  loop  stability  boundary  and  hence  the  disturbance  level 
at  which  the  system  may  be  expected  to  diverge.     From  Fig.    11,    one 
sees  that  a  very  stable  airframe  results  in  an  autopilot  which  cannot 
be  driven  unstable  whereas  an  unstable  airframe  aggravates  the  dis- 
turbed stability  problem  causing  the  autopilot  to  diverge  for  very  small 
inputs.  .         ■        '  .      \   V  .' 
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Figures   13  and  14  illustrate  the  influence  of  actuator  delay  on 
the  stability  boundaries  and  hence  on  the  stability  of  the  overall 
system.     As  can  be  readily  seen,    the  safe  operating  area  decreases 
rapidly  as  T     increases.     Also,    the  frequencies  on  the  acceptable 
operating  region  of  the  boundary  are  much  lower  for  large    T  . 
Lower  frequencies  result  in  larger  chatter  amplitudes,    and  con- 
sequently more  wear  and  tear    on  the  system  hardware.     Small  horse- 
power actuators  with  delays  of  O.Ol   seconds  are  in  existence  but  any 
appreciable  reduction  in  this  figure  will  be  difficult. 

SMALL  SIGNAL  ANALYSIS 

The  steady  state  mode  of  operation  for  the  system  under  analysis 
is  a  high  frequency  limit  cycle.       This  limit  cycle  linearizes  the  non- 
linear element  so  that  a  linear  analysis  may  be  carried  out  by  assign- 
ing a  pseudo-linear  gain  term  to  the  switch.     Such  an  analysis  is  valid 
only  so  long  as  the  limit  cycle  is  not  interrupted  during  any  portion 
of  the  response. 

In  developing  the  stability  surfaces,    the  describing  function  or 
gain  at  the  oscillation  frequency  of  the  nonlinear  element  is  taken 
as     N.       It  can  be  shown  (Z)  that  the  effective  gain  of  the  switch  for 
small  control  signals  is  just  N/2, 

K^=    N/2  ■  (14) 

With  this  simple  relation,    one  may  determine  the  gain 
appropriate  for  a  linear  small  signal  analysis  directly  from  the 
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the  stability  surface.     With  a  pseudo-linear  gain  for  the  switch  thus 
established,    conventional  linear  theory  will  yield  the  response  of 
the  autopilot  for  small  signals. 

TIME   DOMAIN  ANALYSIS 
Switch  Characteristic 


The  preceding  parts  of  this  paper  have  shown  a  method  for  the 
steady  state  analysis  of  the  nonlinear  autopilot  being  discussed. 
Further,    the  small  signal  response  of  the  system  can  be  determined 
from  this  analysis  coupled  with  conventional  linear  theory.       Large 
signals,    however,    interrupt  the  oscillation  and  linear  theory  no  longer 
s\iffices  to  determine  the  system  response. 

With  the  loop  operating  in  a  steady  linnit  cycle,    the  switch  may 
be  characterized  as  sketched  below.     The  smaller  linear  portion  of 
the  characteristic  results  from  the  limit  cycle  linearization  of  the  switch. 
For  nnost  cases,    the  linear  region  is  so  small  relative  to  m,    that  the 
switch  characteristic  may  be  taken  as  that  of  an  ideal,    infinite  gain 
linniter. 
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Optimum  Response 

For  a  bang-bang  or  switch  type  controller,    there  exists  an 
optimum  response  to  an  input  or  disturbance.       Such  a  response 
reduces  the  error  to  zero  in  a  minimum  possible  time.       It  can  be 
shown  that  this  optimum  response  occurs  when  the  two  rules  below 
for  switching  the  limited  variable  are  satisfied. 

1.  The  linnited  variable  is  at  its  maximum 
absolute  value  throughout  the  response, 

2.  The  linriited  variable  changes  sign  (n-  1  ) 
tinnes  where    n    is  the  number  of  poles. 

According  to  Schmidt,    (4),    these  rules  were  derived  for  a  plant  with 
real  poles  but  are  applicable  to  systems  with  complex  poles  provided 
the  damping  is  not  extremely  light. 

From  the  definition  of  optimum  response,    the  error  and    n-  1 

of  its  derivatives  must  be  zero  at    T      ,     the  response  or  solution 

m  ^ 

time.     Note  that  for  step  inputs  the  error  derivatives  are  the  feedback 
quantity  derivatives.       Thus  the  optimum  response  parameters  (the 
reversal  times     Tj,    and  T2,    and  the  response  time     T      )  may  be  found 
by  writing  the  error  and  error  derivative  equations  and  solving  them 
simultaneously.       This  method,    of  course,    depends  on  rules  1  and  2 
to  determine  these  error  equations. 

A  much  more  elegant  and  direct  method  of  finding  the  optimum 
response  parameters  is  to  use  rules  1  and  2  and  the  theory  of  Laplace 
transfornns.     To  this  end,    consider  the  simple  control  loop  shown  in 
the  sketch. 
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Assuming  that  the  input    n      is  a  unit  step  at    t  =  0,    and  that  the 
system  is  designed  to  give  optimum  response,    then  the  error 
e  (t)     obeys  the  conditions: 


e(t)    =    0 
e(t)    =    0 


t    <     0 
t    >    T 


m 


(15a) 
(I5b) 


Now    e(t)     is  a  truncated  time  signal  (i.  e.  ,    zero  for  all  but  a 
finite  time  interval  )  and  must,    according  to  the  theory  of  Laplace 
transforms,   have  a  transform  which  is  an  entire  function.     Using 
rules  I  and  2  and  the  sketch,    one  may  write  the  Laplace  transform 
of  the  system  error.     Forcing  this  function  to  be  entire  yields  the 
optimvim  switching  parameters. 

An  analysis  of  the  type  outlined  above  yields  the  optimum 
switching  or  reversal  times  and  the  optimum  response  times. 
Armed  with  these  data,    it  is  a  simple  matter  to  calculate  the  actual 
response  as  well  as  the  error  and  its  derivatives  during  the  trans- 
ient portion  of  the  response.     Furthermore,    one  may  use  these  data 


to  synthesize  an  optinnum  system. 


Optinaiim  Control  System  for  a  Sinmple  Airframe 
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Derivation  of  Switch  Time  Equations.       To  a  first  approxima- 
tion the  airframe  transfer  function*  is 


n       _    AE 
5      ^       s^ 

Using  this  approximation  of  the  airframe  and  the  methods  outlined, 
an  optimum  system  design  will  be  consumated. 

Since  the  plant  is  third  order,    rule  2  requires  two  switch 
reversals.     From  this  and  rule  I,    the  optimum  form  of  switching 
is  that  shown  in  the  sketch  below.    -^ 
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Using  delay  functions,     5  (t)     is  written  in  terms  of  the  unknown 
switch  times  as  ■  ■ 

6  (t)  =  B  Q  -  2u(t-  Ti)  +  2u(t-  T2)  -  u(t-  T^r]  (l7) 

which  has  Laplace  transform 


6(s)  = 


r        ,  -TjS    ^     ,  -T2S  -T^S~f 

[_l-2e      *      +2e      "-e      "^J 


(18) 


On  the  next  page  is  shown  that  portion  of  the  system  associated 


*    This  representation  has  often  been  used  in  actual  system 
analysis  but  is,    of  course,    of  restricted  accuracy. 
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with  the  analysis  to  follow. 


6(s) 


n,(s) 


^ 


n 


•6(s) 


From  this  sketch  one  writes 


€(s)     =     n^(s)   -  n(s) 


=  n(s)   -^^    6   (s) 

,3 


n 


c 
s 


s" 
BAE 


(19a) 
(19b) 


1  _  Ze-'^^U  2e-'^2^ 


-T     s 
e      m 


J       (19c) 


where  the  latter  form  of  the  equation  assumes  that    n      is  a 

c 

step.     According  to  Schmidt  (4)  the  optimum  solution  is  given 
if  Eq.   (l9c)  is  such  that  there  are  no  poles  in  the  finite  s-plane. 
This  condition  is  satisfied  provided: 

1/3 


n 


Ti    =    ( 


2ABE 
T2    =     3Ti 

T 


) 


n 


nn 


=    4T,  =  4   (- 


2ABE 


) 


1/3 


(20a) 
(20  b) 
(20  c) 


From  these  equations  note  the  following  pertinent  items: 

1.     The  response  time  varies  with    n       such 

c 

that  an  8g  maneuver  takes  twice  as  long 
as  a  Ig  maneuver. 
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2. 

Since    A    and    E    are  each  proportional  to  q, 

the  response  time  varies  as    q          .       Thus 

it  takes  considerably  longer  to  pull  a    2g 

manuever  at  50,  000  feet  altitude  than  it  does 

at  sealevel. 

3. 

The  response  time  decreases  as  the  control 

surface  rate  limit  is  increased.     However, 

actuator  power  consumption  and  chatter  ampli- 

tude both  vary  linearly  with  this  rate  limit. 

4. 

The  switching  is  symmetrical  and  thus  the  con- 

trol surface  rate  is  negative  for  the  same  length 

of  time  as  it  is  positive.     Thus,    of  course,    the 

final  tail  position  is  the  same  as  the  initial  position. 

It  should  also  be  noted  that  although  the  response  time  varies  with    q 

and    n   ,     thi 
c 

s  is  usually  acceptable  from  a  guidance  system  point  of 

view  provided  certain  bounds  are  observed. 

Synthesis  of  Optimum  System.        The  feedback  quantities  in  this 

system,  are 

n,     n,  9  ,     6  ,       These  may  be  written  (at  the  switch  times) 

as: 

n{T 

,.          _     BAE      3    _     ''c                                                             \^.. 
l)                  6       ^'    -      12                                                           ^^^^> 

n(T 

,)          -    ^f   Tf-B^/^    ^1'    n//^                                    (21b) 

e(T.)         =     f     Tf    ~b1/3    q-^/^    n//^                                 (21c) 

e(T,)          =     BE     T,  ~b2/3    ql/3    n^^/^                                         (21d) 
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From  symmetry  of  optimum  switching: 

n(T2)     =     lln^  111  (22a) 

n(T2)     =     A(T,)  (22b) 

0(T2)     =    e(Ti)  (22c) 

e(T2)    =     -e(Ti)  (22d) 


the  end  of  the 

response: 

t    =    T 

m 

n    =    n 
c 

•                      •• 

G    =   e 

=    A   =    5    = 

0 

Expressions  for  the  feedback  quantities,   assuming  an  optimum 

response,    are  given  by  Eqs.   (2])  and    Eqs.   (22).       These  data  must 

now  be  used  to  select    K*    ,   K..  ,    K.,   and  K   ,     such  that  optimum 

b  Q         n  n  *■ 

switching  actually  occurs.     The  optimum  switching  criteria  are 
simply  constraints  on  the  actuating  signal  given  as: 

m(Ti)    =    0  (23a) 

m(T2)    =    0  (23b) 

9m(Ti)     <  Q  3m(T2)      >Q 

9  t  9  t 

The  partial  derivatives  which  insure  proper  direction  of  switching 

were  considered  in  writing  the  equation  for    6(t)  and  hence,    are 

automatically  satisfied.     From  Fig.    15,    the  general  form  of  the 

actuating  signal  is: 

m(t)  =  n^(t)   -  K^  n(t)   -  K^    0(t)   -  Kg   e(t)   -  K^  n(t)  (24) 


32 


Substituting  Eqs.   (Zl)  andEqs.   (22)   into  Eq.   (24)  gives,    respectively: 

fr^^         .  ^^""c  2/^  1/3  2/3 

m(T,)    =    0    =    -^   -  K^K.  R,n^        -  K..R3n^       -  K- R,n^ 

(25a) 


n  2/3  2/3 

e 


m(T,)     =     0    =    ^   .  K^K.  R,n^         -  K-  R,n^         ,  ^.^.n^' 


:  (25b} 

In  these  equations,     K^  =  i  and  the  R's  are  functions  of  A,    B. 
and  E,    and  are  constants  for  a  given  flight  condition.     The 
definition    K^  =  K^K^     is  also  used  in  these  equations. 

Inspection  of  Eqs.   (25)   shows  that    n^    appears  to  various 
powers  in  the  terms  of  the  equations.     Thus  these  equations  can- 
not be  satisfied  for  all  n^    with  constant  values  for  the  feedback 
gains.     If,    instead  of  linear  gains,    one  uses  the  cube  root  of  the 

error  and  the  square  root  of  the  rate  feedback  terms,    the    n 

c 

terms  in  the  equations  may  be  divided  out.     Rewriting  Eqs.  (25) 
in  the  form 

m(T,)  =  0  =^  -  K.    (K^R,+  R3)   n/^'  -  KgR3n/^'  (26a) 

m(T2)  =  0  =^  -  K^    (K^Ri+Rz)  n/    \  K^Rsn/^^  (26b) 

and  using  the  suggested  non-linear  gains  gives 

Iln       ^/^ 
m(T,)  =  0  =  {-^-  )  .  K^^^^1-^R^     n/    '  -  Kg  R3n'/'  (27a) 

m(T3)  =  0  =  (%-)'^lK.    ^n^R;T^3    n/^%K,R3n /^'  (27b) 


sim- 


These  equations  no  longer  involve    n^    and  may  be  solved 
ultaneously  for    K^     and  Kg  .     These  values  would  then  yield 
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a  system  which  is  optimum  for  all    n   .     Note,    however,   that  the 

c 

equations  still  contain  the  aerodynamic  parameter     q    and,    hence, 
the  above  optimum  solution  is  valid  only  at  the  particular     q    for 
which  the  solution  was  taken. 

An  approximation  to  the    q    variation  can  be  made  by  taking 
appropriate  combinations  of  the  feedback  parameter  pair    9     and    n 
and  of  the  pair  6   and  G  .       The  use  of    6   feedback  was  not  considered 
in  the  frequency  domain  analysis  presented  earlier,    so  only  an  in- 
dication of  the  results  in  the  time  domain  will  be  given  here. 

Referring  to  Eq.   (26a),    the  coefficient  of  the    K*      term  is 

K  Ri  +  R2 

r    '■ 

which  in  terms  of    q    becomes 

2/3  -1/3 

K^q  +  q 

where  the  constant  multiplier  has  been  discarded.     It  is  desired  to 

make  this  ternn  as  near     q    independent  as  possible  using  the  simple 

relation  above  (a  polynomial  in    q    could  be  introduced  at  the  expense 

of  great  complication  of  the  hardware).       The  parameter    K      is  selected 

by  eq\iating  the  above  term  for  two  different  values  of    q    and  solving 

for    K   .       Thus 
r 

2/3  -1/3  -1/3  1/3 

Kj.  qi  +  qi  =  qz  +  K^  qz  (28) 

for  a  reasonable  set  of    q,    and    qj    (although  not  necessarily 
optimum)  one  finds    K^  =  1.5,    is  satisfactory.       Using  an  approach 


34 


similar  to  that  above,    one  finds  that  a  suitable  ratio,     K^^  of  5 
feedback  to    G     feedback  is  .0  3. 

Analog  Data  for  Optimum  System.       Actual  response  traces, 
taken  from  an  analog  simulation  of  the  above  system,    are  shown 
in  Figs.    15  and  l6.       For  this  simulation,    diode  function  genera- 
tors were  used  to  approximate  the  nonlinear  functions.     The  analog 
computer  traces,    on  which  is  shown  the  theoretical  optimum  re- 
sponse tinnes,    demonstrate  the  validity  of  the  above  approach. 

The  above  system  has  adequate  response  over  a  large    q 
range  without  the  use  of  air  data  scheduling.     Furthermore,    the 
system  is  faster  than  a  normal  linear  system  designed  around  a  com- 
parable airframe.       Introduction  of  the  nonlinear  gain  terms  precludes 
the  possibility  of  driving  the  system  unstable  with  large  inputs. 

An  important  feature  of  an  optimum  response  is  the  absence 
of  overshoot  without  the  consequent  long  response  time  associated 
with  over  damped  linear  systems.     This  feature  is  evident  in  Figs.   16 
and  17. 

Switch  Time  Equations  for  Neutrally  Stable  Airframe 

The  preceding  discussion  has  considered  a  simple  approximation 
to  the  actual  airframe.     Despite  the  approximation,    switch  time  data 
obtained  from  this  analysis  are  in  surprisingly  good  agreement  with 
the  results  of  more  accurate  analyses. 

Almost  every  nnissile  airframe  has,   at  some  point  in  the  flight 
reginae,    a  condition  of  neutral  stability.       This  simply  implies  that  a 
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perturbation  in  the  angle  of  attack  causes  no  restoring  moment  to 
correct  for  the  disturbance.       From  a  control  point  of  view,    a 
neutrally  stable  airframe  is  desirable  since  such  an  airframe  is  easy 
to  control  and  at  the  same  time  does  not  und\ily  aggravate  the  loop 
statility  problem.     It  is  never  possible,   however,   to  maintain  a  given 
degree  of  airframe  stability  over  the  entire  flight  regime.     Con- 
sequently,   it  is  necessary  that  the  autopilot  loop  be  designed  to  cope 
with  a  range  of  values  of  the  airframe  stability  parameter  C.     The 
two  most  important  cases     C  =  0,     and    C  >  0  are  considered  below. 
A  neutrally  stable  airframe  has  transfer  function: 


n  AE 


6  s^  (s  +  AK  ) 

Since  this  is  a  third  order  plant,     the  equation  for    5(s)  is; 


(29) 


5(s)  =  -f    [l-2e-T^%    2e-T^^     -  e'^-^  ]  ( 30) 

Using  Eqs.   (l8)  and  (l9),    the  error  transform  can  be  written  as: 

€(s)       !^    -     BAE  [l  -  2e-^^^+  Ze""^^^     -  e-^^s]  ^^^^ 

^  s^(s  +  AK) 

This  function  has  no  poles  in  the  finite  s-plane  if: 

0    =     2Ti  -  2T2+  T^  (32a) 

AKTi  A'KTz  AKT 

0=1-  2e  +2e  -e  ^  (32b) 

n  K  T     ^ 

_c _     rp  2  2  m  ,        , 

BE    "     ^2    -  Ti —  (32c) 

Useful  literal  expressions  for    Tj,    Tz,   and  T        cannot  be 

m 
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obtained  from  these  equations  necessitating  a  numerical  machine 
method  of  solution. 

Switch  Time  Equations  for  Unstable  Airframe 

An  unstable  airframe  is  defined  by    C>  0  and  may  be  caused 
by  center  of  gravity  shifts  with  changing  fuel  load,    center  of  pres- 
sure shifts  or  any  number  of  other  causes.       It  was  previously  shown 
that  an  unstable  airframe  results  in  a  control  loop  which  diverges  for 
much  smaller  inputs  or  distrubances.      A  rather  surprising  result  of 
the  analysis  to  follow  is  that  the  response  time  for  the  loop  with 
C  >  0  is  actually  longer  than  for  the  C  =  0  case. 

With  C  >  0     the  airframe  transfer  is 


where 


n 


AE 


s(s^  +  AKs  -  C) 


^      AK-f  '/(AK)^+  4C 


AE 


az    = 


AK   -    ^fiAK)^+  4C 


s(s  +  ai)(s  +  az) 


a;,    >   0 


a,  <    0 


The  error  transform  for  this  case  is 


n 


-BAE 


e(s)  =-£ _ 

^  As+  ai)(s  +  az) 


(33) 

{34a) 
(34b) 


•T,s 


■T,s 


-  I  -  2e  +  2e 


•T     s 
m 


(35) 
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For  optimum  switching,    one  gets: 

aiTi  aiT2         ajT 

0  =  1  -  2e  +  2e  -  e 

0  =   1  -  2e  +  2e  -  e 


&-    2  T,  -  2  T.  -  T^ 


m 


m 


(36a) 
(36b) 
(36c) 


Solution  of  these  equations  was  also  carried  out  using  machine 
computation. 

Optimum  Switch  Time  Data 

Optimum  switch  time  data  for  the  flight  conditions  in  Table  2 
are  presented  in  Figs.    18  through  21. 


Table  2.     Aerodynamic  parameters  used  in  finding 
optimum  switch  time  data. 


Flight 

Control  Surface 

condition 

A 

D 

E 

K 

Rate 

Limit  B 

1 

1.0 

0 

229.0 

1.  1 

200 

8 

1.0 

49.0 

229.0 

1.  1 

200 

4 

1.0 

0 

229.0 

0 

200 

5 

0.  116 

0 

26.0 

1.7 

200 

7 

0.116 

5.53 

26.0 

1.7 

200 

6 

0.116 

0 

26.0 

0 

200 

These  figures  illustrate  the  small  differences  in  switch  times  associated 
with  the  various  airframes.       As  would  be  expected,    the  unstable  air- 
frame requires  the  first  switch  point  to  occur  early,    otherwise,    the 
angle  of  attack  woxild  get  uncontrollably  large.       The  rather  unexpected 
result,    however,    is  the  longer  response  time  associated  with  the 
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unstable  airframe. 

Synthesis  Procedure 

With  the  simple  airframe  considered  earlier,    it  was  fairly- 
easy  to  write  literal  expressions  for  the  switching  times  and  switch- 
ing equations.       The  more  exact  airframe  representation,    however, 
requires  the  use  of  numerical  miethods  as  outlined  below. 

As  before,    the  switching  criteria  are 

m(Ti)    =    0  (37a) 

m(T2)    =  0  (37b) 

9  m(Ti)  9    m(Tz) 

—5-1 "'  5~t  " 

provided  there  is  no  delay  time.       To  achieve  optimvim  response, 
the  limited  variable    6    must  actually  change  sign  at    Tj    and    T2. 
If  the  actuator  has  delay,    then  the  actuating  signal  must  anticipate 
the  reversals  by    t    seconds.       Hence,    the  switching  criteria  in  the 
presence  of  delay  become* 

m(Ti  -  t)  =  0  (38a) 

m(T2  -  t)  =  0  (38b) 

9    MTi  -  t)       ^   Q    d   m(Tg  -  t)     ^q 

at  at 

The  partial  derivatives  are  satisfied  by  the  expression  for  6(s) 
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and  the  other  constraints  may  be  written  as: 

m(Ti-  t)  =  0  =  n^-  K^n(Ti-  x)   -  K^K^  A(T,-  t)   -  K^e(Ti-  t)  -Kge(Ti-T) 

(39a) 

m(T2-T)  =  0  =  n^-  K^n(T2-  t)   -  K^K^  ACT^-  t)  -  K- e(T2-  t)  -  Kg0(T2-T) 

,  -.■,..  (39b} 

The  feedback  quantities  may  easily  be  foiond  from  the  equations: 

n(t)  =  L-^  p(s)   G(s)3  (40a) 

A(t)  =    -^  n(t)  (40b) 

6(t)  =  n(t)  K  +  ^^  (40c) 

e(t)  =  E  6(t)+  C   ^^  (40d) 

6  (t)  =  B  |j  -  (t  -  Ti)  u  (t  -  Ti)  +  (t  -  Tz)  u  (t  -  T2)] 

With    K      and    K       specified,   Eqs.   (39)  maybe  solved  simultan- 
r  n 

eously  for  the  optim\im  values  of    K-      and    K.  .       To  this  end,   a 
digital  computer  program  was  prepared  to  yield: 

1.  The  optimum  switch  times  for  C=  0  and  C  >  0. 

2.  The  solutions  of  Eqs.  (40)  and  (4l)  at  t  =  Tj-  t  and  t  =  T2  -  t. 

3.  The  simultaneous  solution  of  Eqs.  (39)  for  optinnum  values 

of    Ki     and    K"    • 
0  0 

4.  The  computation  and  plotting  of    n(t)  from  Eq.  (40a). 

Figure  22  shows  the  values  of    K*     ,   as  determined  fronn  Eqs. 
(39),    plotted  versus  the  command  acceleration    n   .       Data  are  shown 
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for  the  two  extreme  flight  conditions.       Notice  that  the  optimum    K; 
values  for  the  low    q    condition  are  very  nearly  equal  to  those  for 
the  high    q    case.       This  is  due  directly  to  the  use  of    n    feedback 
(which  introduces  the  parameter    K    )  in  combination  with    9     feed- 
back to  render    Kt     independent  of    q.       The  results  shown  in  Fig.    21 
indicate  that  1.5  is  a  good  choice  for    K    . 

The  optimum  gain  data  shown  in  Figs.    22  and  23  are  strong 

fvinctions  of    n   .       To  remove  this  dependence,   at  least  two  nonlinear 
c 

gain  functions  must  be  introduced.       To  find  appropriate  functions  for 
these  gains,    one  must  apply  numerical  methods  to  the  data  obtained 
from  Eqs.   (39)  and  Eqs.   (40)  .       This  operation  was  not  carried  out 
since  the  data  obtained  are  valid  only  for  a  specific  set  of  airframe 
parameters  and  are  thus  of  little  general  use.       Results  obtained 
analytically  for  the  simple  airframe  indicate  the  general  nature  of 
these  nonlinear  gains.       In  selecting  nonlinear  gains,    it  should  be 
noted  that  the  small  signal  values  are  determined  by  the  frequency 
domain  stability  analysis. 
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NOMENCLATURE 

5 

control  surface  rate 

B 

control  surface  rate  limit 

m(t) 

actuation  signal 

G(s) 

forward  path  transfer  function 

H(s) 

feedback  path  transfer  function 

q  = 

1/2  kpM^  dynamic  pressure 

N 

describing  function  of  non-linear  element 

a 

sine  wave  amplitude  at  input  to  non- linear  element 

CO 

o 

angular  oscillation  frequency 

U(Ni,w^ 

real  part  of  characteristic  equation 

V(Npa)J 

imaginary  part  of  characteristic  equation 

K 

n 

acceleration  feedback  gain 

K. 
n 

acceleration  rate  feedback  gain 

'^6 

body  angular  rate  feedback  gain 

Kg 

body  angular  acceleration  feedback  gain 

K 

r 

''■J'^i 

A 

g's  per  degree  angle  of  attack 

D 

g's  per  degree  control  surface  incidence 

C 

body  angular  acceleration  per  degree  angle  of  attack 

E 

body  angular  acceleration  per  degree  control  surface 

incidence 

K 

angle  of  attack  rate  due  to  lateral  acceleration 

T 

actuator  delay  time 

A 

AE  -  DC 
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NOMENCLATURE    concl. 

00        =    -=-         accelerometer  natural  frequency 
na         T  ^ 

a 

w        =    -=r  gyro  natural  frequency 

"^  g 

C  accelerometer  damping  ratio 

C  gyro  damping  ratio 

n  connmand  acceleration 

c 

5  control  surface  incidence 

Wz  high  frequency  oscillation  frequency 

t*3 1  low  frequency  oscillation  frequency 

Kj-  gain  of  the  non-linear  element  for  small  signals 

e  system  error 

Tj  first  reversal  time  of  optimum  response 

Tz  second  reversal  time  of  optimum  response 

T  optinnum  response  time 

a  =  T^  T^ 
a      g 

^  =  24    T    T2+  2;    T    T^ 
a    a   g  g    g      a 

7=  T^  +  T^  +  44    T    4    T 
a  g  a    a    g    g 

^  =  24    T     +24    T 
a    a  g    g 

0    =  2KA4    T    +    E 
a    a 

TT=  KAT^  +  24    T    E 
a  a    a 

4)  =  ET^ 
^  a 

H^   =  2  4     T 
g     g 


_               ""I 

' 
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ctf 

<^H 

fe 

■P 

•P 

b 

12 
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5 

^ 
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^ 
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Fig.  7  K©  stability  boundaries  for  various  K^ 
at  high  q  condition. 
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Fig.  8  K^  stability  boxmdaries  for  various  K^ 
at  low  q  condition. 
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Fig.  9    Kq  stability  boundaries  for  various  Kq 
at  high  q  flight  condition. 
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Fig,  10    Kg> stability  boundaries  for  various  K^ 
at  a  lov  q  flight  condition. 
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^ 


airframe  stability  on  system  stability. 
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Fig,  12  I^  stability  boundary  illustrating  influence 
of  airframe  stability  on  system  stability. 
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APPENDIX  A 
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The  missile  airframe,   which  is  an  important  part  of  the  system 
discussed  in  this  paper,    is  of  the  tail  control  type.       Smiall  control 
surfaces,    located  aft  of  the  center  of  gravity,    are  used  to  generate 
turning  nnoments  on  the  missile  body.     This,    in  turn,    rotates  the 
body  to  an  angle  of  attack  such  that  the  body  forces  provide  the  re- 
quired lateral  acceleration. 

A  representation  of  this  airfranne  for  control  system  analysis 
purposes,    due  to  Amsler,    is  shown  in  the  simple  block  diagram  of 
Fig.    I-A.       The  representation  used  here  assumes  that  perturbations 
about  a  trim  condition  are  defined  by  the  linear  equations: 

Q     =    Eb+  C  a 

a     -    9-7  =  9-  nK 

n    =    Aq-  +  D5 

For  certain  trim  conditions  and  very  large  perturbations,    the  airframe 
behavior  may  depart  significantly  from  linearity.       A  study  of  system 
behavior  for  nonlinear  aerodynannics  is  beyond  the  scope  of  this  paper. 
Working  directly  from  Fig.    I-A,    one  finds  the  following  airframe 
transfer  functions: 

n  Ds^  +    A 


6  s(s2+  AKs  -  C) 

9       _  Es+  KA 

6  s(s2  + AKs  -  C) 

A  =  AE  -  DC 


(A-1) 


(A- 2) 
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In  general,     D    is  very  small  and  may  be  neglected.       For  various 

airframe  stability  conditions  and  for    D=  0    Eq.   (A-l)  becomes: 

-                                                                   ^  Ota  Die  airirame/ 

5            s(s2+  AKs  -  C) 

'  ,                                                                         ^iNeutraiiy  stable  airirame/ 

5            s^s  +  AK) 

"^  -                    ^^                   ^'  >  °         fTT     1   n       ■    r          ^ 

^Unstable  airframe} 

5          s(s  +  ai)(s  +  slz)        a2  <  0 

A    = 

C^      qs/mg 

D    = 

C         qs/mg 
^^6 

E    = 

C          qsd57.3/I 

6 

C    = 

C          qsd57.3/I 

a 

K    = 

57.3xg/V^ 

^M 

=    moment  coefficient 

^N    = 

=    normal  force  coefficient 

^C^r 

""n^ 

N 
9  a. 

a 

9C^, 

\ 

N 
-       36 

V 

9C,, 

""^^ 

M 
da 

a 

9C^. 

^=-5 

M 
'6b 
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q    =    dynamic  pressure 

s    =    reference  area 

d    =     reference  length 

m    =    airframe  mass 

I  s:   moment  of  insrtia  about  center  of  gravity 

V        =    missile  velocity 
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APPENDIX  B 
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Stability  b  oundary  data  presented  in  this  paper  were  computed 
using  the  IBM  1620  computer  operating  with  a  Fortran  program. 
The  source  statments  for  this  problem  along  with  computer  variables 
definition  in  terms  of  system  variables  are  listed  below. 

1  READ  2,  A,  B,  C,  E,  FKBAR,  TG,  DG,  TA,  DA,  FKl 

2  FORMAT  (F7.  3,  FT.  3,  F7.  0,  F7.  0,  FT.  1,F7.5,F7.  1,F7.5, 

F7.  1,    F7.  l) 

3  READ  4,FK2,FK3,TPR,FK 

4  FORMAT  (F7.5,F7.5,F7.5,F7.5  ) 

5  PUNCH  6 

6  FORMAT  (47H     A      B      C      E      KBAR      K2      K3) 

7  PUNCH  8,   A,  B,  C,  E,  FKBAR,  FK2,  FK3 

8  FORMAT(F7,3k  F7.3.F7.  3»F7.3,   F7.  3,  F7.  3,  F7.  3) 

V  =    TA*  *  2 
X    =    TG**2 

Y  =    DA*TA 

Z    =    DG*TG 
ALPHA    =    V*X 

BETA    =     2.*(Z*  V+Y*X) 

GAMMA    =     V+X+4.  ♦  Y*  Z) 

RPIO    =     2.*(Y+Z) 

DELTA    =    A*E-B*C 

THETA    =    2.  *  FKBAR*  DELTA*  Y+E 

PI    =    FKBAR*  DELTA*  V+ 2.  *Y*E 

PHI    =    E*V 

FMU    =!    2.*Z 

FK4    =    FK*FK2 

K    =    0 

W    =     1. 
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30     WW    =     W*  *  2 

TMPA    =     X*WW-1. 

TMPB    =    B*  WW-DELTA 

TMPC    =     PHI*WW-THETA 

TMPD    =    PI*  WW-FKBAR*DELTA 

TMPE    =     TMPA*  TMPB 

TMPF    =    FMU*W*TMPB 

TSIN    =    SINF(W*TPR) 

TCOS    =     COSF(W*  TPR) 

All    =     -(BETA+A*  FKBAR*  ALPPIA)*  WW*  *  3+(FHCH-A*  GAMMA* 
FKBAR-C*  BETA)  *  WW*  WW 

A12    =     -(A*FKBAR-RHO*C)*WW 

Al  =  All  +  A12 

Bl    =    TMPC*  WW*  TCOS- W*  TMPD*  TSIN 

Cll  =  FK4*  W{ TMPF*  TCOS+ TMPE*  TSIN)  -FK2* 
(TMPD*  TCOS+W*  TMPC*  TSIN  ) 

C12  =  FKl*  (TMPE*  TCOS-TMPF*  TSIN) 
CI   ^     Cll  +  C12 

A22  =  W*( -ALPHA*  WW*  *3+(GAMMA+A*FKBAR*  BETA- C 
*  ALPHA  )  *  WW*  WW  ) 

A23  =  -( l.+RHO*  A*  FKBAR-C*  GAMMA)*  WW*  W-C*W 
A2  =  A22  +  A23 

B2  =  W*  (-TMPD*  TCOS-W*  TMPC*  TSIN) 

C22  =  FK4*  W*(TMPE*  TCOS- TMPE*  TSIN)  +  FK2*  W* 
(-TMPC)  *TCOS+  TMPD*  TSIN) 

C23  =  -FKl*  (TMPF*  TCOS+TMPE*  TSIN  ) 

C2  =  C22  +  C23 

DNOM  =  A.1*  B2-A2*  Bl 

IF(DNOM)   100,    200,100 

100    RECN=  (B1*C2-B2*C1  )/dNOM 

IF(RECN)  200,200,  110 
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110     FK5  =  (A2*  C1-A1*C2)  /DNOM 

B3  =  FK3*Bl+FKl*(TMPE*TCOS-TMPF*TSIN) 

C3  =  FK*  W*(TMPF*TCOS+TMPE*TSIN)-(TMPD*TCOS 
+  W*TMPC*TSIN) 

B4  -  FK3*B2-FKl*(TMPF*TCOS+TMPE*TSIN) 

C4  =  FK*  W*(TMPE*TCOS-TMPF*TSIN)+(W* 
( -  TMPC  )  *  T  COSf  TMPD*  TSIN  ) 

NDOMl  =  Al*  C4-A2*  C3 

IF(FK5)   200,  200,  120 

120  IF(DNOMl)    121,200,121 

121  RECNl  =  (B4*C3-B3*  C4  )/    DNOMl 

122  FK6  =  (A2*  B3-A1*  B4)  /  DNOMl 

123  K  =   1 

PUNCH  94,  W.RECN,  FK5,RECN1,FK6 
94    FORMAT  (F6.  1,  F12.  5,  F12.  5,  F12.  5,  F12.  5  ) 
150    IF(W-14.5)    152,152,153 

152  W  =  W  +  .  5 
GO    TO    30 

153  IF(W  -   1.65/TPR)   154,300,300 

154  W  =  W  =  5. 
GO    TO    30 

200    IF(k)   300,  150,  300 
300    CONTINUE 

GO    TO    3 

END 
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Computei'  variable 

System  variable 

TA 

^a 

TG 

\ 

DA 

^a.     ■ 

DG 

'^:      .] 

FKBAR 

K 

FK4 

K. 

n 

FK 

K 

r 

FK2 

% 

W 

u 

B 

D 

TPR 

T 

FK3 

•% 

FK5 

Computed    K" 
0 

FK6 

Comiputed    K* 
0 

> 
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This  thesis  presents  a  method  of  analyzing  a  nonlinear  control 
system  to  determine  both  stability  and  response.     The  stability 
analysis  is  carried  out  in  the  frequency  domain,    and  the  response 
analysis  is  treated  in  the  time  domain.       For  the  specific  system 
of  the  type  discussed  in  this  paper,    these  nnethods  are  complemen- 
tary.      Typical  numerical  data,    obtained  by  using  these  methods, 
are  presented. 

The  autopilot  on  which  the  analysis  is  carried  out  consists  of 
an  aerodynannic  vehicle  with  control  surfaces  driven  by  bang-bang 
type  actuators  with  transportation  lag.       The  aerodynamic  loop  is 
closed  with  two  sensing  instruments  in  each  steering  plane;    an 
angular  rate  gyro  and  a  lateral  accelerometer.       These  instruments 
are  coupled  with  appropriate  mixing  and  shaping  networks  to  yield  a 
system  which  does  not  require  air  data  scheduling. 

Describing  function  methods  are  applied  to  the  autopilot  loop  to 
yield  an  approxinaate  characteristic  equation.       From  this  equation, 
stability  surfaces  in  a  parameter  three  space  which  define  the  pos- 
sible oscillatory  modes  are  derived.      Root  locus  methods  are  used 
to  interpret  these  surfaces  and  to  determine  what  parts  of  the  three 
space  are  acceptable  operating  regions.       It  is  shown  that  the  upper 
portion  of  the  stability  surface  is  an  unsafe  boundary  and  must  be 
avoided.       The  lower  portion  of  the  surface  defines  a  steady  limit 
cycle  oscillation  which  linearizes  the  control  loop  for   small  signals. 
Using  the  coordinates  of  an  operating  point  on  this  part  of  the  surface 


and  linear  theory,    one  may  find  the  small  signal  response  of  the 
autopilot  loop. 

Numerical  stability  data  for  a  set  of  typical  flight  conditions 
and  feedback  parameters  are  shown.       In  addition,    sets  of  data  are 
presented  to  show  the  influence  of  airframe  stability  and  actuator 
delay  on  the  overall  loop  stability. 

For  large  signals  the  limit  cycle  oscillation  is  suppressed 
throughout  most  of  the  transient  period.       During  this  time,    the  non- 
linear element  behaves  as  a  switch  or  limiter  and  is  treated  as  such. 
Optimum  response  for  a  switch  type  controller  is  assumed  to  exist 
and  the  theory  of  Laplace  transforms  used  to  find  this  optimum 
solution.       Using  these  data,    feedback  networks  may  be  designed 
which  will  yield  the  assximed  optimuni  response.       For  a  simplified 
airframe  representation  analytical  expressions  are  developed  for 
the  feedback  parameters.      A  numerical  procedure  for  designing  the 
feedback  networks  for  a  complex  airframe  is  suggested.       The  in- 
fluence of  airframe  stability  on  the  optimum  bang-bang  solution  is 
pointed  out. 


